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Abstract 

Based on the low-energy effective Hamiltonian with the generalized factor- 
ization, we calculate the new physics contributions to branching ratios and 
CP-violating asymmetries of the charmless hadronic decays B s — * h\hi in 
the standard model and the general two-Higgs doublet models (models I, II, 
and III). Within the considered paramter space, we find the following, (a) In 
models I and II, the new physics corrections are always small in size and will 
be masked by other larger known theoretical uncertainties, (b) In model III, 
the new physics corrections to the branching ratios of those QCD penguin- 
dominated decays B s — > K°T](\ K + K~* ,etc, are large in size and insensi- 
tive to the variations of M H + and N% . For tree- or electroweak penguin- 
dominated decay modes, however, the new physics corrections are very small 
in size, (c) For B s — » K + K~* and other seven decay modes, the branching 
ratios are at the level of (1 — 3) x 10~ 5 and will be measurable at the future 
hadron colliders with large b production, (d) Among the studied thirty nine 
B s meson decay modes, seven of them can have a CP-violating asymmetry 
Acp larger than 20% in magnitude. The new physics corrections are small 
or moderate in magnitude, (e) Because of its large and N£ stable branching 
ratio and CP violating asymmetry, the decay B s — ► K + K~* seems to be the 
"best" channel to find CP violation of B a system through studies of two-body 
charmless decays of B s meson. 
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I. INTRODUCTION 



In B experiments, new physics beyond the standard model (SM) may manifest itself, for 
example, in following two ways 0,0]: (a) decays which are expected to be rare in the SM are 
found to have large branching ratios; (b) CP-violating asymmetries which are expected to 
vanish or be very small in the SM are found to be significantly large or with a very different 
pattern with what predicted in the SM. These potential deviations may be induced by the 
virtual effects of new physics through loop diagrams. 

The observation of many two-body charmless hadronic B u d meson decays by CLEO, 
BaBar and Belle 0-0], the successful start of the asymmetric B factories at SLAC and KEK, 
and the expectation for large number of events of B u ^ meson decays to be accumulated at 
B factories and other hadron colliders stimulated the intensive investigation for various B 
decay channels. The two-body charmless hadronic decays B u ^ — ► /ti/12 [ where h\ and h 2 
are the light pseudo-scalar (P) and/or vector(V) mesons ] have been studied, for example, 
in Refs. @-|nfl. 

It is well known that the low energy effective Hamiltonian is the basic tool to calculate the 
branching ratios and Acp of B meson decays. The short-distance QCD corrected Lagrangian 
at NLO level is available now ] 14, 15] , but we do not know how to calculate hadronic matrix 



element from first principles. One conventionally resort to the factorization approximation 
16| . However, we also know that non-factorizable contribution really exists and can not be 



neglected numerically for most hadronic B decay channels. To remedy the naive factorization 
hypothesis, some authors [|TT|,p~0|,p7T|] introduced a phenomenological parameter AT? (i.e. the 



effective number of color) to model the non-factorizable contribution to hadronic matrix 
element, which is commonly called the generalized factorization. Very recently, Cheng et al. 
[]T8|| studied and resolved the controversies on the gauge dependence and infrared singularity 
of the effective Wilson coefficients C e J^ |19|] by using the perturbative QCD factorization 
theorem. 

Unlike the B u ^ meson, the heavier B s meson can not be produced in CESR, KEKB and 
PEP-II e + e~ colliders. Only upper limits on decay rates of several charmless hadronic B s 
decays are current available from LEP collaborations |2"0,2"T], such as B s — > K + K~, K + tt~, 
7c i] and B s — > 7777, while most of them are far beyond the theoretical predictions. However, 
it is expected that many B s decays can be seen at the future hadron colliders with large b 
production. Recent theoretical studies and experimental measurements about the mixing of 
B® — B® can be found in Refs. [22,23|. The early studies of two-dody charmless hadronic 
decays of B s meson can be found in Refs. p^ , p5| . Based on the framework of generalized 
factorization, Tseng analyzed the exclusive charmless B s decays involving r/'\ while 



Chen, Cheng and Tseng |TI| calculated the branching ratios of thirty nine charmless two- 
body decays of B s meson. It is found that the branching ratios of 777/ - 1 and several other 
decay modes can be as large as 10~ 5 and measurable at future experiments. 



In a recent work ||27|| , we made a systematic study for the new physics contributions to 
the branching ratios of seventy six B u ^ — > h\h 2 decay channels in the framework of the 
general two-Higgs-doublet models (2HDM's). In this paper we extend the work to the case 
of B s meson. In additional to the branching ratios, we here also calculate the new physics 
contributions to the CP- violating asymmetries Aqp of charmless hadronic decays B s — > h\h 2 
induced by the new gluonic and eletroweak charged-Higgs penguin diagrams in the general 
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2HDM's ( models I, II and III). Using the effective Hamiltonian with improved generalized 
factorization [|TI| , we evaluate analytically all new strong and electroweak penguin diagrams 
induced by exchanges of charged Higgs bosons in the quark level processes b — ■> qV* with 
q G {d, s} and V G {gluon, 7, Z}, and then combine the new physics contributions with their 
SM counterparts and finally calculate the branching ratios and CP-violating asymmetries 
for all thirty nine exclusive B s — > h\h2 decay modes. 

This paper is organized as follows. In Sec. |IJ, we describe the basic structures of the 
2HDM's and examine the allowed parameter space of the general 2HDM's from currently 
available data. In Sec. [TTT], we evaluate analytically the new penguin diagrams and find the 
effective Wilson coefficients with the inclusion of new physics contributions, and present 
the formulae needed to calculate the branching ratios B(B — > hih 2 ). In Sec. [V] and [V], 
we calculate and show numerical results of branching ratios and CP-violating asymmetries 
for thirty nine B s decay modes, respectively. We focus on those decay modes with large 
branching ratios and large CP-violating asymmetries. The conclusions and discussions are 
included in the final section. 



II. THE GENERAL 2HDM'S AND EXPERIMENTAL CONSTRAINTS 



The simplest extension of the SM is the so-called two-Higgs-doublet models |28[ . In such 
models, the tree level flavor changing neutral currents(FCNC's) are absent if one introduces 
an discrete symmetry to constrain the 2HDM scalar potential and Yukawa Lagrangian. Lets 
consider a Yukawa Lagrangian of the form 



£y = rEQwfaUjjt + r}%Qi,L<f>iD jiR + &Qi,i.<hUj,ii + &Qi,L<hDj, R + h.c 



* = 1,2) 
) with % 



are the two Higgs doublets of a two-Higgs-doublet model, 0^2 = 2T 2 0i 2 , 
[1,2,3) are the left-handed isodoublet quarks (right-handed up-type 
quarks), Dj R are the right-handed isosinglet down- type quarks, while rj^f and ^j- = 
1,2,3 are family index ) are generally the non-diagonal matrices of the Yukawa coupling. 



n, 
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TUi, one 



By imposing the discrete symmetry: X 
obtains the so called model I and model II. 

During past years, models I and II have been studied extensively in literature and tested 
experimentally, and the model II has been very popular since it is the building block of 
the minimal supersymmetric standard model. In this paper, we focus on the third type of 
the two-Higgs-doublet model 



usually known as the model III p9| , |30| . In model III, no 
discrete symmetry is imposed and both up- and down-type quarks then may have diagonal 
and/or flavor changing couplings with 0i and 02- As described in | 29[] , one can choose a 
suitable basis (H°, H 1 , H 2 , i7 ± ) to express two Higgs doublets. The H ± are the physical 
charged Higgs boson, H° and h° are the physical CP-even neutral Higgs boson and the A 
is the physical CP-odd neutral Higgs boson. After the rotation of quark fields, the Yukawa 
Lagrangian of quarks are of the form 



£¥ 1 = 'iJh.lO.rj.H + vSQi,L<PlDj,R + ^Qu.O-J'j.H + &Qi,L<hDj,R + H.c. 



(2) 



where r^- D correspond to the diagonal mass matrices of up- and down-type quarks, while 



the neutral and charged flavor changing couplings will be |p9 
on the &' D couplings as the Ref. |29 



We make the same ansatz 
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(3) 



^1,2,3) are 



where Vckm is the Cabibbo-Kobayashi-Maskawa mixing matrix 
the generation index. The coupling constants Ay are free parameters to be determined by 
experiments, and they may also be complex. 

In model II and setting 1 < tan/? = v^/vi < 50 favored by experimental measurements 
POfl , the constraint on the mass of charged Higgs boson due to CLEO data of b — ► 57 is 



> 200 GeV at the NLO level [32 



For model I, however, the limit can be much 
weaker due to the possible destructive interference with the SM amplitude. For model III, 
the situation is not as clear as model II because there are more free parameters here p9| , |33l . 
In a recent paper Chao et al. studied the decay b — ► 57 by assuming that only the 
couplings Aft = \\tt\e t6t and A&& = \\bb\e t6b are non-zero. They found that the constraint 
on Mh+ imposed by the CLEO data of b — > 57 can be greatly relaxed by considering the 
phase effects of X u and A^. From the studies of Refs. p4| , |35| , we know that for model III 
the parameter space 



Kj '- 
|A«| 



0, for ij 7^ tt, or bb, 
= 0.3, |A»| = 35, 9 = (0 C 



30 c 



M H+ = (200 ± 100)GeV, 



(4) 



are allowed by the available data, where 9 = 6^ — 6 tt . 

From the CERN e + e _ collider (LEP) and the Fermilab Tevatron searches for charged 
Higgs bosons p6| , the new combined constraint in the (M^+tan/5) plane has been given, 
for example, in Ref. ||20|| : the direct lower limit is M H + > 77 GeV, while 0.5 < tan/3 < 60 
for a relatively light charged Higgs boson with ~ 100 GeV. Combining the direct and 
indirect limits together, we here conservatively consider the range of lOOGeV < M#+ < 300 
GeV, while take Mh+ = 200 GeV as the typical value for models I, II, and III. For models 
I and II we consider the range of 1 < tan (3 < 50, while take tan f3 = 2 as the typical value. 

III. EFFECTIVE HAMILTON! AN IN THE SM AND 2HDM'S 



The standard theoretical frame to calculate the inclusive three-body decays b — > sqqf\ is 



based on the effective Hamiltonian [|T5inrTin~3 



Gp I 



H eff (AB = 1) = -J rECj (v ub V: s Q] + V^Qf) - V tb V* 
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^ CjQj + CgQg 

3=3 



(5) 



Here the first ten operators Q\ — Q10 can be found for example in Refs. [n,l3", 27[, while the 
chromo-magnetic operator reads: 



Q 9 = i^m b sa^(l + ls)T a af3 bpG% 



8tt 2 



(6) 



For b — > dqq decays, one simply make the replacement s — > d. 
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where a and (3 are the 577(3) color indices, T% ( a = 1, ...,8) are the Gell-Mann matrices. 
Following Ref. [fHql , we do not consider the effect of the weak annihilation and exchange 
diagrams. 

The coefficients Cj in Eq.(^j) are the well-known Wilson coefficient. Within the SM 
and at scale Mw, the Wilson coefficients Cx(Mw), ■ ■ ■ , C\q{Mw) and C g (Mw) have been 



given for example in Refs. [fL4] , |15(| . By using QCD renormalization group equations, it is 



straightforward to run Wilson coefficients Ci(Mw) from the scale /i = 0(M\y) down to the 
lower scale n = 0(mb). Working consistently to the NLO precision, the Wilson coefficients 
Ci for i = 1, . . . , 10 are needed in NLO precision, while it is sufficient to use the leading 
logarithmic value for C g . 



A. New strong and electroweak penguins 



For the charmless hadronic decays of B meson under consideration, the new physics 
will manifest itself by modifying the corresponding Inami-Lim functions Cq(x), D (x), E (x) 
and E' (x) which determine the coefficients Cs(Mw), ■ ■ ■ , C\q{Mw) and C g (Mw)- These 
modifications, in turn, will change the SM predictions of the branching ratios and CP- 
violating asymmetries for decays B s — > h\h>2 under study. 

The new strong and electroweak penguin diagrams can be obtained from the correspond- 
ing penguin diagrams in the SM by replacing the internal W lines with the charged-Higgs 



H + lines. In Ref. ||27|| , we calculated analytically the new Z -, 7- and gluon-penguin dia- 



grams induced by the exchanges of charged-Higgs boson H + , and found the new Co, D , E , 
and E'q functions which describe the new physics contributions to the Wilson coefficients 
through the new penguin diagrams, 
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J(y t )\\tt\ 2 -K(y t )\\ tt \ bb \e 



w 



with 



H(y) 

m 

J(y) 

K(y) 



+ — rr— My] 



72(1 -yf 
IQy - 29y 2 + 7y 3 

36(1 - yf 
2y + 5y 2 - y 3 

4(1 -yf 
-3y + y 2 



+ 
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3y 2 
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logM, 
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(7) 

(8) 
(9) 
(10) 

(11) 
(12) 
(13) 
(14) 



4(1 -y)* 2(1 -y)3 

where x t = m 2 /M 2 v , y t = m 2 /Mj I+ , and the small terms proportional to m 2 /m 2 have been 
neglected. In models I and II, one can find the corresponding functions Cq, D ,E and E' Q 
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by evaluating the new strong and electroweak penguins in the same way as that in model 
III: 



r 1 — r 11 — Xt 
G °" G ° "Stan 2 /? 
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(15) 
(16) 
(17) 
(18) 
(19) 



where y t = mf/Mf I+ , tan/5 = Vifv\ where v\ and Vi are the vacuum expectation values of 
the Higgs doublet 4>i an d 02 as defined before. 

Combining the SM part and the new physics part together, the NLO Wilson coefficients 
Ci(M w ) and C g (M w ) can be written as 



C X {M W ) 
C 2 (M W ) 
C 3 (M W ) 

C 4 (M W ) 

C 5 (M W ) 

C 6 (M W ) 

C 7 (M W ) 
C 8 (M W ) 
C 9 (M W ) 
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sin 2 fl 



C g (M } 



w 



\(E'{x t ) + Er 



(20) 
(21) 

(22) 

(23) 

(24) 

(25) 

(26) 
(27) 

(28) 
(29) 



where x t = mf/M^, the functions B (x), C (x), D (x), E (x) and E' Q are the familiar 
Inami-Lim functions |37] in the SM and can be found easily, for example, in Refs. [|I~4"| , |38" . 
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Since the heavy new particles appeared in the 2HDM's have been integrated out at the 
scale Myy, the QCD running of the Wilson coefficients Cj(M^) down to the scale \x = 0{m^) 
after including the new physics contributions will be the same as in the SM: 



C{ i i) = U(ji,M w )C(M w ), 
C g ^)=r, 1A ^C g {M w ) + Y^hri a ' 



(30) 
(31) 



where C(M W ) = (Ci(M w ), . . .,C 10 (M W )) T , U(fi,M w ) is the five-flavor 10 x 10 evolution 
matrix at NLO level as defined in Ref. rj = a s (Mw)/oi s ([i), and the constants hi and 
Oj can also be found in Ref. [ |I~4]j . 

In the NDR scheme and for SU (3)c, the effective Wilson coefficients f] can be written as 





cf ff 



1 + ^ K + 7y log- 



a 



a. 



C j + 7^-A' i (C t + C p + C g ) + -^B> i Q 



ij 



24tt 



8tt 



(32) 



where A[ = (0, 0, -1, 3, -1, 3, 0, 0, 0, 0) T , B[ = (0, 0, 0, 0, 0, 0, 1, 0, 1, 0) T , the matrices f v and 
7v contain the process-independent contributions from the vertex diagrams. The matrix 7^ 
and fy have been given explicitly, for example, in Eq.(2.17) and (2.18) of Ref. Note 
that the correct value of the element (f NDR ) 66 and {f ndr)88 should be 17 instead of 1 as 



pointed in Ref. [[39 



The function C t , C p , and C g describe the contributions arising from the penguin diagrams 
of the current-current Qi,2 and the QCD operators Q3-Q6, and the tree-level diagram of the 
magnetic dipole operator QsG; respectively. We here also follow the procedure of Ref. [T^ 
to include the contribution of magnetic gluon penguin. The functions C t , C p , and C g are 
given in the NDR scheme by HTT| , p~3 



C* 



C n 



c P 



a, 



G(m q ) - G{m h ) 



C 3 + 



C u 
10 
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i=u,d,s,c,b 



(C 4 + C 6 ), 



I + ^G(m u ) + ^G(m c ) 



(C 1 + 3C 2 ), 



(33) 

(34) 

(35) 
(36) 



with X q i = V q 'bV*, q , and C|" = C g {ji) +C5. The function G(m) can be found, for example, in 
Refs. [|l^,[27]]. For the two-body exclusive B meson decays any information on k 2 is lost in the 



factorization assumption, one usually use the "physical" range for k 2 |TTHl3 
Following Refs. [ITT1HT3] we take k 2 = ml/2 in the numerical calculation. 



< 1.2 < "<K 
2 



2 In the improved generalized factorization approach [Q, these effective coefficients are renormal- 
ization scale- and scheme-independent, gauge invariant and infrared safe. 
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B. Decay amplitudes in the BSW model 



Following Ref. the possible effects of final state interaction (FSI) and the contri- 



butions from annihilation channels will be neglected although they may play a significant 
rule for some decay modes. The new physics effects on the B decays under study will be 
included by using the modified effective coefficients di (i — 3, . . . , 10) as given in the second 
entries of Table | and Table || for the model III. The effective coefficients a, in models I and 
II are not shown explicitly in Table | and Table |T[ In the numerical calculations the input 
parameters as given in Appendix and Eq.(|4|) will be used implicitly. 

With the factorization ansatz |Tj|, the three-hadron matrix elements or the decay am- 
plitude < XY\H e ff\B s > can be factorized into a sum of products of two current matrix 
elements < X\ Jf|0 > and < Y\J 2fJ ,\B s > ( or < Y\ Jf |0 > and < X\J 2fJl \B s >). The explicit 
expressions of matrix elements can be found, for example, in Refs. JTBLBDI . 

In the B rest frame, the branching ratios of two-body B meson decays can be written as 

B{B S -> XY) = r Bs ^M_\M(B S - XY)\ 2 (37) 
for B s — > PP decays, and 



8nM Bs 



3 



B(B S -> XY) = r Bs ^]r 72 \M{B s -> XY)/(e ■ p B )\ 2 (38) 



8vrM- 



for B s — > PV decays. Here t{B®) = lA93ps P0 |, p B is the four-momentum of the B meson, 
My and e is the mass and polarization vector of the produced light vector meson respectively, 
and \p\ is the magnitude of momentum of particle X and Y in the B rest frame, 



H = 2M B V [M * " {Mx + My)2][M2b ~ (Mx ~ My)2] (39) 

For B s — > VV decays, the situation is more involved. One needs to evaluate the helicity 
matrix element H\ =< Vi(\)V2(\)\H e ff\B) > with A = 0, ±1. The branching ratio of the 
decay B — > V1V2 is given in terms of H\ by 



B(B S - VM) = r Bs (|^ |2 + \ H+ll 2 + 1^2) (4Q) 



M 

MM 

where \p\ has been given in Eq.(|39D. The three independent helicity amplitudes H , H + i 
and H_i can be expressed by three invariant amplitudes a, b, c defined by the decomposition 



10 a ^ 



a9 ^ + MW/^ + MM e ^ PlP 



(41) 



where pip and M\p are the four momentum and masses of Vj.,2, respectively, p = p\ + p 2 is 
the four-momentum of B meson, and 

H±i = a± cVx 2 - 1, H = -ax - b (x 2 - l) (42) 

Ml-M 2 -Ml , t , 

x = — ? i ^ 43 

2. /;. /j V 7 
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For individual decay mode, the coefficients a, b and c can be determined by comparing the 
helicity amplitude H\ =< Vi(\)V2(\)\H e ff\B a > with the expression fl4i"l) . 

In the generalized factorization approach, the effective Wilson coefficients Cf will ap- 
pear in the decay amplitudes in the combinations: 

a 2i _ 1 ^C 2 f_ 1 + ^, a* = C$ + ^fg, (* = 1,..., 5) (44) 

where the effective number of colors N^ s is treated as a free parameter varying in the range 
of 2 < Aff < oo, in order to model the non-factorizable contribution to the hadronic 
matrix elements. Although A^ ff can in principle vary from channel to channel, but in the 
energetic two-body hadronic B meson decays, it is expected to be process insensitive as 
supported by the data (Ej. As argued in ref. |T7[, Nf(LL) induced by the (V - A)(V - A) 
operators can be rather different from N° S (LR) generated by (V — A)(V + A) operators. 
Since we here focus on the calculation of new physics effects on the studied B meson decays 
induced by the new penguin diagrams in the two-Higgs doublet models, we will simply 
assume that N^ S (LL) = N° S (LR) = N^ s and consider the variation of N° s in the range of 
2 < A^f < oo. For more details about the cases of N^ S (LL) ^ N° S (LR), one can see for 



example Ref. |12||. We here will not consider the possible effects of final state interaction 
(FSI) and the contributions from annihilation channels although they may play a significant 
rule for some decay modes. 

Using the input parameters as given in Appendix, and assuming k 2 = ml/ 2, M H + = 200 
GeV, the theoretical predictions of effective coefficients a$ are calculated and displayed in 
Table | and Table |lj for the transitions b^d(b^d) and b — > s (b — > s), respectively. 
For coefficients a 3 , . . . , a w , the first and second entries in tables (0,0) refer to the values of 
ctj in the SM and model III, respectively. 

Compared with Ref. JT^], the effective coefficients given here have two new features: 

• The effective Wilson coefficients here are not only renormalization scale- and 
scheme-independent, but also gauge invariant and infrared safe. 

• The contribution due to the chromo- magnetic dipole operator Q„ has been included 



here through the function C g as given in Eq . fl36|) . For the penguin dominated decay 
channels, operator C g will play an important role. 

• The coefficient a\ and a 2 remain unchanged in 2HDM's since the new physics consid- 
ered here does not contribute through tree diagrams. 

• The new physics contributions are significant to coefficient a 4 and a^, but negligibly 
small to coefficients 03,5 and 07-10- 

All branching ratios here are the averages of the branching ratios of B and anti-5 decays. 
The ratio 5B describes the new physics correction on the decay ratio and is defined as 

5B(B - XY) - B{Bs XY)NP - B{Bs XY)SM U5) 
5B(B S - XY) B{Bs - xy)SM (45) 
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IV. BRANCHING RATIOS OF B s MESON DECAYS 



Using the formulas and input parameters as given in last section and in Appendix, it 
is straightforward to find the branching ratios for the thirty nine B s — > PP, PV, VV decay 
channels. In the numerical calculations, we use the decay amplitudes as given in Appendix 
A, B and C of Ref. [12| directly without further discussions about details. 



Following Ref. the hadronic charmless B meson decays can be classified into 

six classes: the first and last three classes correspond to the tree-dominated and penguin- 
dominated amplitudes, respectively. 

• Class-I and Class-II decays are dominated by the external and internal W-emission 
tree-diagrams, respectively. Examples are B s — »■ K + ix~ , K Q ix , 

• Class-Ill: the decays involving both external and internal W-emissions. But this class 
does not exist for the B s decays. 

• Class-IV and Class-V decay modes are governed by effective coefficients 04,6,8,10 an d 
a 3,5,7,9, respectively. Examples are B s — > K + K~ , irrf \ ■ ■ ■ . 

• Class- VI decays involve the interference of class-IV and class-V decays. 

In tables pTT[ - |VT| , we present the numerical results of the branching ratios for the thirty nine 
B s — > PP, PV, VV decays in the framework of the SM and models I, II and III. Theoretical 
predictions are made by using the central values of input parameters as given in Eq.(|J) 
and Appendix, and assuming A = 0.804, A = 0.22, p = 0.16, 77 = 0.34, M H + = 200GeV, 
9 = 0°,30°, tan/3 = 2 and N* s = 2, 3, 00 in the generalized factorization approach. The 
fc 2 -dependence of the branching ratios is small in the range of k 2 = ml/2 ± 2 GeV 2 and 
hence the numerical results are given by fixing k 2 = ml/2. 

The SM predictions for all B s decay modes as listed in tables [TTT| and [IV] are agree well 
with those given in Ref. ]12| . The effect of changing ry and including the new contribution 
from the chromo- magnetic operator Q g in the SM is not significant. 

For decay modes involving B s — > K* or B s — > transitions, we use two different 



set of form factors: the Bauer, Stech, and Wirbel (BSW) form factor [16[ and the light- 
cone sum rule(LCSR) form factor as given explicitly in Appendix. For decay modes 
B s — > 7r°0, (f)T]' , p°(j), ujcf) and B s — > <fxft, the variation of the branching ratios induced by 
using different set of form factors is about a factor of 2, but small or moderate for all other 
decay modes. 

From numerical results, we see the following general features of new physics corrections: 

• In model III, the new physics corrections to QCD-penguin dominated decay modes, 
such as B s — > K°7]^\ rf'rp\ K°K°, etc., are large in size and insensitive to the vari- 
ations of the mass M H + and N* s : from 30% to 130% w.r.t the SM predictions for 
both cases of 9 = 0°, 30°. For tree-dominated or electroweak penguin dominated decay 
modes, however, the new physics corrections are very small in size: SB < 5%. 

• In models I and II, the new physics corrections to all B s — > hih 2 decay modes are 
always small in size within the considered parameter space: less than 10% and 20% 
in model I and II respectively, as shown in tables and VI. So small corrections will 
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be masked by other larger known theoretical uncertainties. Variation of tan/3 in the 
range of 2 < tan (3 < 50 can not change this feature. 

In model III, the new gluonic penguins will contribute effectively through the mixing 
of chromo- magnetic operator Q g with QCD penguin operators Q 3 — Q 6 , as shown 
in Eq . (|32|) . The C e J^ will strongly dominate the new physics contributions to B s 
meson decays. The branching ratios for all thirty nine decay modes have a very weak 
dependence on 9 in the range of 0° < 9 < 30°. 



As pointed in Refs. fl2y41|l , the decays 

B a ^r}n,7}'n,r}p,7}'p,^m,(f>p (46) 

do not receive any QCD penguin contributions, and are predominately governed by ag and 
hence N£ s insensitive. In 2HDM's, this remain to be true because the new physics corrections 
to coefficients a7_io are negligibly small as shown in tables | and [TI], and therefore, the new 
physics contributions to these decay modes are also very small: < 2%. As suggested in Ref. 
|T2|j , a measurement of these six decay modes can be utilized to fix the parameter ag. It is 



clear that the inclusion of new physics contributions in the 2HDM's does not change this 
picture. 

For the decays 

B s -> our], 4*1^ , K<f>, K*<f>, 4>n, (47) 

the SM electroweak penguin corrections are in general as important as QCD penguin effects, 
and very sensitive on . The new physics corrections to these decay modes in the model III 
also have a strong dependence on the variation of N° s : 5B = —20% — 110% for 2 < N* s < oo. 
As illustrated in Fig. p], for example, the branching ratio of B s — > <pr] decay has a moderate 
M#+-dependence, but a strong N* s dependence. For Fig. |l](a) and 0(b), we set N° s = 3 and 
Mh+ = 200 GeV, respectively. The four curves correspond to the theoretical predictions in 
the SM (dotted curve), model II (dot-dashed curve), model III with 9 = 0° (solid curve) and 
9 = 30° (short-dashed curve), respectively. 

Among the thirty nine charmless two-body hadronic B s decays, we find that only 7 (8) 
of them have branching ratios at the level of 10~ 5 in the SM (model III): 

B s -> K + K~, K°K°, TjTj, rjlrf\ K + p, K + *p~, 00. (48) 

Among these eight decay modes, the new physics correction to the Class-I decay mode 
B s — > K + p~ and K + *p~ are very small, from — 2%to 1%. For remaining six decay modes, 
the new physics enhancement are significant: from ~ 50% to ~ 130% and insensitive to 
the variation of N% . These decay modes will be measurable at the future hadron colliders 
with large b production P2| . In Figs. and |, we plot the mass and N* s dependence of the 



branching ratios of B s — > K + K~~ and r\vj decay modes. 

After inclusion of new physics contributions in the models I, II and III, the patterns 
observed in Ref. [0 remain unchanged, 

T(B S -> K+K-) > T(B 8 -> K + K*-) > T(B 8 -> K* + K*-) > T(B S -> K + *R-), 
T(B S -> K°K°) > T(B S -> K°K*°) > T(B S -> K*°K*°) > F(B S -> K*°K°). (49) 
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Recently, large decay rates for B+ — > K + r] f and Bd — > K°rj' decays have been reported 
by CLEO and BaBar collaborations The CLEO measurement of B® — > K°r]' decay is 

B(B^ — > K°rf) = (89ltil ± 9) x 10~ 6 , which is larger than the branching ratios of B — > Kn 
decays by a factor of 3 to 5. For B s decays, the decay modes B s — > rjrj' and B s rj'rj' are 
the analogue of Bd — > K°rj' decay, and are expected to have large branching ratios. From 
Table [IH[ one can see that the SM predictions of the branching ratios B(B S — ► rjrj') and 
B(B S — > rj'rj') are indeed large, but in comparable size with other six decay modes listed 
in Eq.(^). The new physics enhancement to these two decay modes are significant in size, 
~ 70% in model III, as illustrated in Fig. [| After inclusion of new physics contributions, 
we find numerically that 

B(B S -> rjrj') sa (23 - 33) x 10~ 6 , (50) 
B(B S -> rj'rj') « (12 - 16) x 10~ 6 . (51) 

These theoretical predictions will be tested by the future experimental measurements. 

For decays B — > K + K~* and B — > K + *K~*, they have relatively large decay rates and 
weak Mh+ and N* dependence. In Figs. |] and we plot the mass and N* s dependence 
of the branching ratios B(B S —>■ K + K~*) and B(K + *K~*). It is easy to see that the new 
physics contributions in the model III to these two Class-IV decays are significant ( ~ 70% 
) in size and insensitive to the variations of and N% s . 



V. CP- VIOLATING ASYMMETRIES OF B s MESON DECAYS 



In Ref. [25], Du et al. studied the branching ratios and CP- violating asymmetries for 



decay modes B s — ► K~tt + , K + K~, K°tc , cfxfi and K°<j). Recently, Ali et al. [52] estimated the 
CP-violating asymmetries in seventy six charmless hadronic decays of B u and Bd mesons. 
The calculation of the CP- violating asymmetry Acp for B s meson decays are theoretically 
very similar with those of the Bd meson decays. For more details about the theoretical 



aspects of CP- violating asymmetries in B u> d h\h,2 decays, one can see Ref. and 
reference therein. In this section, we calculate the CP- violating asymmetries of B s — > hih,2 
decays in the framework of the SM and the general two-Higgs-doublet models. We focus 
on evaluating the new physics effects on Acp for thirty nine B s decay channels induced by 
charged Higgs penguin diagrams appeared in the general two-Higgs-doublet models. 

In models I and II, one does not expect sizable changes in Acp of B s decays since 
there is no any new phase introduced when compared with the SM. In model III, although 
the introduce of a new phase 6 played an important role in relaxing the constraint on the 
parameter space of model III due to the CLEO measurement of B — > X s ^ decay as studied 
in Ref. p4| , we still do not expect dramatic changes for the pattern of the CP- violating 



asymmetries of B s decays under consideration because this phase may alter the theoretical 
prediction of Aqp through loop diagrams only. 

Analogous to the Bd meson decays, the time dependent CP asymmetry for the decays of 
states that were tagged as pure B® or B® at production is defined as 

Acp{t) -T(B0( t )^f) + T(B0(t)^f) {b2) 
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Following Ref. ||42|| , the neutral B® (B®) decays can be classified into three classes ac- 



cording to the properties of the final states / and /: 

• Class-1: B® — > /, B® — > /, and the final states / or / is not a common final state 
of B® and B®, for example, B° s — > K + ir~. The CP- violating asymmetry for class-1 
decays will be independent of time: 

T(B s ^f)-r(B° s ^f) 
•A-cp = =n — (53) 

r(B.° - f) + r(Bj - /) 

in terms of partial decay widths. 

• Class-2 and 3: — > (/ = /) with / CP = ±/ (class-2) or / CP ^ ±/ (class-3), the 
time-integrated CP asymmetries are of the form 

^ p -l + ^l + |A C pp- 2 l + a;n + |A C pP (54) 



with 



^V to </lg e// lB°> 
P " ^</|^//|50> (55) 



where x = AM B o/T B o « 20 is the preferred value in the SM [25| for the case of 
B® — Bg mixing |. Contrary to the B^ meson decay where x ~ 0.73, it is easy to see 
that the parameter x for B° s decays is very large. The first and second term in Eq.flSip 
is strongly suppressed by 1/x 2 and 1/x, respectively. We therefore do not expect large 
CP- violating asymmetries Aqp for the class-2 and class-3 B Q S decays. This expectation 
is confirmed by the numerical results given below. 



In Tables |VI1| and - |VIII| , we present numerical results of CP- violating asymmetries Aqp 



for thirty nine B s — ► h\hi decay channels in the SM and 2HDM's, using the input parameters 
as given in Appendix, and assuming that k 2 = ml/2, p = 0.16, r\ = 0.34, M H + = 200GeV, 
6 = 0°, 30°, and Nf* = 2, 3, oo. We show the numerical results for the case of using BSW 
form factors only since the differences induced by using the BSW or LCSR form factors are 
small for almost all B s decay modes. 

Among thirty nine B s decay modes studied, we find that seven of them have CP- violating 
asymmetries larger than 20% in the SM and model III: 

B s -> K°*n°, K°p°, K°*r) { '\ K + K~*, K°*p°, K°*u . (56) 

All these seven decay modes are belong to the CP-class-1 decay modes. On the other 
hand, all twenty four class-2 and 3 decay modes have small CP-violating asymmetries only, 
\Acp\ 5%, mainly due to the strong suppression of 1/x 2 as shown in Eq. 



From Ref. M, the upper limit is x = AM B o/T B o > 15.7 at 95%C.L. 
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In models I and II, the new physics corrections on Acp for almost all B s decay modes 
studied here are negligibly small as can be seen from Table VIII and Figures In model 



III, the new physics correction is varying from channel to channel, as illustrated in Table 
VI I| and Figures 



• For B s — > K + K decay, the new physics correction to its Acp is very small in size 
and insensitive to the variations of and 9. 

• For B s — > K + K~* decay, the new physics correction to its Acp is moderate in size, 
from —20 to —40% with 0° < 9 < 30°, and insensitive to the variations of N* . 

• For B s — > K°*T]' and three remaining decays given in Eq.(^), the size and the sign of 
the new physics corrections strongly depend on both N% and 9. 

• For B s — > 777/, (fxf) and several other CP-class-2 and 3 decays, the new physics correc- 
tions can be as large as a factor of 30, but have a very strong dependence on N^ s and 
9. Despite the large new physics correction to these decay modes, their Acp are still 
smaller than five percent because of strong suppression of 1/x 2 . 

For the QCD-penguin dominated B s — > K + K~* decay, its decay amplitude is propor- 
tional to the combination of large and N^ G stable coefficient ai and a 4 |12[ . 



M(B S -> K + K-*) oc [V ub V* ai - V tb V*(a 4 + a 10 )] . (57) 

The imaginary part of M. for b — > s and b — > s transitions are very different, which in turn 
leads to a large Acp- Numerical result indeed shows that this decay has a large and N^ s 
stable CP-violating asymmetry, 

A C p(b\ -> K ± K^*) w -30% (58) 

for 2 < < 00. Another advantage of this decay mode is the large (~ 70%) new physics 
enhancement to its branching ratio B(B S — ► K + K~*) in model III, as illustrated in Fig.[7]. 
Taking into account above facts, this decay mode B s — > K + K~* seems to be the "best" 
channel to find CP violation of B s system through studies of two-body charmless decays of 
B s meson. 

Since the tree-dominated B s — > K + tt~ decay mode has a moderate CP-violating asym- 
metry ( ~ 10% ), a large branching ratio (~ 7 x 10~ 6 ), negligible new physics correction, 
large detection efficiency and a very weak Nfr dependence, we therefore classify this decay 
mode as one of the promising decay channels for discovering the CP violation in B s system. 

For the decay B s — > K°*1]', although the SM prediction of its Acp can be large, but it is 
varying in the range of —60% to 60% due to the strong dependence on , as illustrated in 
Fig||. Another disadvantage for this decay is its small branching ratio, (0.02 — 0.16) x 10~ 6 , 
almost two orders smaller than that of B s — > K + n~ and K + K~* decays. 



4 In general, the detect efficiency for the two-body B meson decays with charged final states is 
larger than that with neutral final states by a factor 2 or 3. 
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For the remaining five decay modes as given in Eq. fl56|) , although the size of their Acp can 
also be as large as 20 to 30%, but these decays can not be "good" channels for discovering the 
CP violation in B s system because of the strong N^r dependence and very small branching 
ratios. 

In Figures |]-(7|, we show the mass and Nfr dependence of Acp for B s — > K + K~ and 
K + K~* decays. In these figures, the dotted and dot-dashed curve refers to the theoretical 
prediction in the SM and model II, while the solid and short-dashed curve corresponds to the 
prediction in the model III for 6 = 0° and 30° respectively. As can be seen from Fig.[7|, the 
CP-violating asymmetry of B s —>■ K + K~* decay are large in size and has weak or moderate 
dependence on M H +, Nf 1 and 6. 



VI. SUMMARY AND DISCUSSIONS 

In this paper, we calculated the branching ratios and CP-violating asymmetries of two- 
body charmless hadronic decays of B s meson in the standard model and the general two- 
Higgs-doublet models (models I, II, and III) by employing the NLO effective Hamiltonian 
with generalized factorization. In Sec.|T|, we defined the effective Wilson coefficients G\' 
with the inclusion of new physics contributions, and presented the formulas needed to cal- 
culate the branching ratios B(B S — ► h\h2). 



In Sec. |TV|, we calculated the branching ratios for thirty nine B s — > h\h,2 decays in the 
SM and models I, II, and III, presented the numerical results in Tables |T^-(VJ and displayed 
the M H + and N^r dependence of several interesting decay modes in Figures |lf|| From the 
numerical results, one can see the following 

• In models I and II, the new physics corrections to the decay rates of all B s — > h\}\2 decay 
modes are small and will be masked by other larger known theoretical uncertainties. 

• In model III, the new physics corrections to QCD penguin-dominated decays B s — > 

are large in size, from 30% to 130% w.r.t the SM predictions, 
and insensitive to the variations of the mass M#+ and N^ s . For tree- or electroweak 



penguin-dominated decay modes as listed in Eq. (f46l) , however, the new physics cor- 
rections are very small in size: 5B < 5%. 

For the decays B s — > t]r] f and B s — > rj'r]', analogue of B^ — > K°r]' decay, the branching 
ratios are large but in comparable size with other six decay modes listed in Eq.(pE5|). 
The new physics enhancements to B(B S — > 7777') and B(B S — > rfr}') are significant in 
size, ~ 70% in model III. 

For decay modes B s 7r°0, 4>tj', p°0, uxp and B s — > 00, the variation of the branching 
ratios induced by using the BSW or LCSR form factors is about a factor of 2, but 
small or moderate for all other decay modes. This feature remain basically unchanged 
after inclusion of new physics contributions. 

For B s — > K + K~ and other decay modes as listed in Eq.(^), the branching ratios 
are at the level of (1 — 3) x 10~ 5 in the SM and model III. These decay modes will be 
measurable at the future hadron colliders with large b production. 
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In Sec. 0, we calculated the CP-violating asymmetries Acp for thirty nine B s — > hihz 



decays in the SM and 2HDM's, presented the numerical results in Tables |VII| - |VIII| and 
displayed the M H + and Nf 1 dependence of Aqp for several typical decay modes in Figures 
|)-|[ From those tables and figures, the following conclusions can be drawn: 

• For almost all B s decay modes, the new physics corrections on Acp are negligibly small 
in models I and II. In model III, the new physics correction is varying from channel 
to channel, and has a strong dependence on the parameter N* s and the new phase 9 
for most decay modes. 

• For twenty four CP-class-2 and 3 B s meson decay modes, their CP-violating asymme- 
tries are small, \Acp\ < 5%, due to the strong 1/x 2 suppression. 

• Among the studied thirty nine B s meson decay modes, seven of them can have a 
CP-violating asymmetry larger than 20% in magnitude. 

• The B s —>■ K + K~* decay has a large and N^ s - and ^-stable CP-violating asymmetry, 
~ —30%, and a large branching ratio. This mode seems to be the "best" channel to 
find CP violation of B s system through studies of two-body charmless decays of B s 
meson. The tree-dominated B s — > K + tt~ decay is also a promising decay channel for 
discovering the CP violation in B s system. 
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APPENDIX A: INPUT PARAMETERS AND FORM FACTORS 



In this appendix we present relevant input parameters. The input parameters are similar 



with those used in Ref. [ |12| . 

• The coupling constants, B meson masses, light meson masses, • • -, are as follows (all 



masses in unit of GeV ) 12.20 



a em = 1/128., a s (M z ) = 0.118, sin 2 6 W = 0.23, G F = 1.16639 x 10~ 5 (Gel/)~ 
M z = 91.188, M w = 80.42, m B o = 5.369, m n ± = 0.140, 
77M, = 0.135, m n = 0.547, m„> = 0.958, m n = 0.770, m w = 0.782, 
= L.Uiy. ;;/ /x= = u. 19 I. in K << = U. UJN. ////, = = 0.892, 



1.019, m K ± = 0.494, m K o = 0.498, m K *± 
0.896, t(B° s ) = 1.493ps, 



(Al) 
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For the elements of CKM matrix, we use Wolfenstein parametrization and fix the 
parameters A, A, p, 77 to their central values: 



A = 0.804, A = 0.22, p = 0.16, 77 = 0.34 



(A2) 



Following Refs. |26|,|T3|, the current quark masses evaluated at the scale p = m b will 
be used in the numerical calculations, 



m b (m b ) = AMGeV, m c (m b ) = 0.95GeV,m s (m b ) = 0A05GeV, 
m d (m b ) = 6AMeV, m u {m b ) = 3.2MeV. 

For the mass of heavy top quark we also use m t = rnt(m t ) = 168GeV. 



(A3) 



For the decay constants of light mesons, the following values are used in the numerical 
calculations (in the units of MeV): 

U = 133, f K = 160, f K . = 221, f p = 210, f u = 195, U = 237, 

fv~fv~ fv' = fv' = fv = — 2 -4> frf = — 6.3, 

= -112, = 137, . (A4) 

where f^ {/) and f^ (l) have been defined in the two- angle-mixing formalism with 6q = 
-9.2° and 8* = -21.2° 131. 



In the calculation we use the following BSW form factors F(0)(in the units of GEV) 



(A5) 



(A6) 



F B_Mr (0) = 0.33, F B -**(0) = 0.274, F^ v (0) = -0.212, F BW (0) = 0.218, 
Kh^i ) = °- 273 ' = 0.236, Af^ K '(0) = 0.232, 

V B -><f>(0) = 0.319, V B ^ K *(0) = 0.2817. 

We use the monopole fc 2 -dependence for form factors, 

fi{k ) — - -pr-/ 2' 
1 — k /m% 

where m* is the pole mass given in fi~6| : 

{m(0"), m(l~), m(l+), m(0 + )} = {5.2789, 5.3248, 5.37, 5.73}, (A7) 

for ub and db currents. And 

{m(0"), m(l~), m(l + ), m(0 + )} = {5.3693, 5.41, 5.82, 5.89}, (A8) 

for sb currents. 

For the decays involving B s — > K* and B s — > <p transitions, we also consider the case 



of using LCSR form factors with the fc 2 -dependence as defined in Ref. [§4 

/(°) 



f(k 2 



l-a(kyMl) + b(kz/Miy 



(A9) 



where the values of /(0) and coefficients a and b have been given in Ref. [fy 
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TABLES 



TABLE I. Numerical values of a« for the transitions b —>■ d \b —>■ d ]. The first, second and third 
entries for 03, ... , aio refer to the values of a, L in the SM and models II and III, respectively. All 
entries for 03, ... , a±o should be multiplied with 10 -4 . 







Nf 


= 2 




Nf 


= 3 


Nf 


= 00 


ai 




0.985 


0.985] 




1.046 


1.046] 


1.169 


1.169] 


02 




0.216 


0.216] 




0.021 


0.021] 


-0.369 


-0.369] 


03 


-10.4 - 


- 19.lt 


-11.5 - 25.7*] 




66.1 


66.1] 


219 + 38.1*' 


221 + 51.4*] 




-33.0 - 


- 19.1* 


-34.1 - 25.7*] 




66.2 


66.2] 


265 + 38.1*' 


267 + 51.4*] 


04 


-349 - 


- 95.3i 


[-354 - 129*] 


-386 


- 102* 


-392 - 137*] 


-459- 114* 


-466 - 154*] 




-463 - 


- 95.3i 


[-469 - 129*] 


-507 


- 102* 


-513 - 137*] 


-596- 114* 


-602 - 154*] 


a 5 


-163 - 


- 19.lt 


-164 - 25.7*] 




-61.5 


-61.5] 


142 + 38.1*' 


144 + 51.4*] 




-186 - 


- 19.lt 


-187-25.7*] 




-61.4 


-61.4] 


187 + 38.1*' 


189 + 51*.4] 




rno 
— OOO " 


nr ' 

- yo.dt 


r CC A A 1 OOnl 

[—044 — LAMl\ 


— 562 


- lUz* 


— 568 — lo (i\ 


— 6(jy — 114* 


—616 — 154*J 




-652 - 


- 95.3i 


[-657 - 129*] 


-683 


- 102* 


-689 - 137*] 


-746- 114* 


-752 - 154*] 


a 7 


5.2 


— 2.5* 


5.1-3.1*] 


4.1 


- 2.5* 


4.0-3.1*] 


2.1 - 2.5*' 


2.0-3.1*] 




5.4 


- 2.5* 


5.3 - 3.1*] 


4.3 


- 2.5* 


4.2 - 3.1*] 


2.2 - 2.5*' 


2.1 - 3.1*] 


a 8 


7.2 


- 1.3* 


7.2 - 1.6*] 


6.9 


- 0.8* 


6.8- 1.0*] 


6.2 


6.2] 




7.4 


- 1.3* 


7.3 - 1.6*] 


7.0 


- 0.8* 


7.0- 1.0*] 


6.3 


6.3] 


a 9 


-85.8 


-2.5* 


-85.9-3.1*] 


-91.7 


- 2.5* 


-91.8-3.1*] 


-103-2.5* 


-104-3.1*] 




-86.4 


-2.5* 


-86.5-3.1*] 


-92.3 


- 2.5* 


-92.4-3.1*] 


-104.1 - 2.5* 


[-104-3.1*'] 


aw 


-16.5 


- 1.3* 


-16.6 - 1.6*] 


0.7 


- 0.8* 


0.7- 1.0*] 


35.2 


35.2] 




-16.6 


- 1.3* 


-16.7- 1.6*] 


0.7 


- 0.8* 


0.7- 1.0*] 


35.4 


35.4] 
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TABLE II. Same as Table | but for b — > s [b — ► s ] transitions. 







N cS 


= 2 




N cS 

1 *c 


= 3 


N cS 


= oo 


(1 1 




0.985 


0.985] 




1.046 


1.046] 


1.169 


[1.169] 


do 




0.216 


0.216] 




0.021 


0.021] 


-0.369 


[-0.369] 


n q 


— 10.9 - 


- 21. 7* 


[_q S — 22 1*1 




66.1 


66.1] 


220 4- 43 3? 


1218 + 44 3*1 




—33.6 - 


- 22. 7* 


—32 5 — 22 2*1 




66.2 


66.2] 


266 + 43 3* 


[264 + 44 3*1 


0, l 


—352 - 


- 108* 


_346 _ iiijl 


—389 - 


- 116* 


_383 — 118*1 


—462 — 130? 

T:L/ l«JU(i 


[_455 _ 133 ? 1 




—467 - 


- 108* 


_4f,n _ 111,1 


—510 - 


- 116* 


_504 — IISi'l 


_5qq _ 13Q ? - 


[_5Q2 — 1337*1 




— 1 64 - 


- 22.7i 


_i fi9 — 22 2*1 

.1.1. i— — . ^ 1 \ 




—61 5 


— 61 51 


1 43 -1- 43 3? 


fl 4D + 44 3*1 

l^ru tit: .Oil 




— 186 - 


- 21.7* 


— 185 — 22 2*1 




—61.4 


—61.4] 


188-1- 43 3? 


[186 + 44 3*1 


u '0 


—541 - 


- 108z 


—535 - 111*1 


-565 - 


- 116* 


_55Q _ 118*1 


—612 — 130* 


[—606 - 133*1 




—654 - 


- 108z 


_64q _ iiijl 


-686 - 


- 116* 


—680 — 118*1 


— 7AQ — 130, 


f— 742 — 133*1 


u f 


5.1 


- 2.8* 


5.2 - 2.8*1 


4.1 


- 2.8* 


4.2 - 2.8*] 


2.0 - 2.8* 


[2.1 - 2.8*] 




5.3 


- 2.8* 


5.4 - 2.8*'] 


4.3 


- 2.8* 


4.4 - 2.8*] 


2.2 - 2.8* 


[2.3 - 2.8*'] 


«8 


7.2 


- 1.4* 


7.2 - 1.4*] 


6.9 


- 0.9* 


6.9 - 0.9*] 


6.2 


[6.2] 




7.4 


- 1.4* 


7.4- 1.4*] 


7.0 


- 0.9* 


7.0 - 0.9*] 


6.3 


[6.3] 


«9 


-85.9 


- 2.8i 


-85.8 - 2.8*'] 


-91.7 


- 2.8* 


-91.6 - 2.8*] 


-104-2.8* 


[-103 - 2.8*] 




-86.5 


- 2.8i 


-86.4 - 2.8*'] 


-92.4 


- 2.8* 


-92.3 - 2.8*] 


-104-2.8* 


[-104-2.8*] 


«lo 


-16.6 


- 1.4* 


-16.5 - 1.4*'] 


0.7 


- 0.9* 


0.7-0.9*] 


35.2 


[35.2] 




-16.7 


- 1.4* 


-16.6 - 1.4*'] 


0.7 


- 0.9* 


0.7 - 0.9*] 


35.4 


[35.4] 



21 



TABLE III. B(B S -> h u h 2 ) (in units of 10~ 6 ) in the SM and model III by using the BSW form 
factors, and assuming k 2 = ml/2, p = 0.16, r] = 0.34, M H + = 200GeV, 9 = 0°, and Nf = 2, 3, oo. 



SM: B Model III: B and 5B[%] 



Channel 


Class 


2 


3 


00 


2 


3 


00 


2 


3 


00 


B s 


-> K + it- 


I 


6.33 


7.14 


8.89 


6.52 


7.35 


9.16 


3.1 


3.1 


3.0 


B s 


T^O 


II 


0.19 


0.08 


0.56 


0.24 


0.14 


0.64 


23.8 


67.1 


14.4 


Bs 


-> K [) 7] 


VI 


0.34 


0.31 


0.79 


0.47 


0.46 


1.00 


38.3 


49.9 


26.8 


B s 


-> K U Tj' 


VI 


0.57 


0.51 


0.77 


0.88 


0.84 


1.17 


53.0 


65.6 


52.3 


B s 


-> K + K~ 


IV 


10.7 


11.7 


14.0 


16.7 


18.5 


22.3 


56.5 


57.6 


59.4 


B s 


n 

— ► 7r u r/ 


V 


0.04 


0.06 


0.11 


0.04 


0.06 


0.11 


1.9 


1.8 


1.3 


B s 


-> 7r u r/' 


V 


0.04 


0.06 


0.10 


0.04 


0.06 


0.11 


1.9 


1.8 


1.3 


B s 


— > 7777' 


VI 


13.8 


15.9 


20.5 


22.5 


25.9 


33.4 


63.8 


63.3 


62.6 


B s 




VI 


6.79 


7.51 


9.08 


11.6 


12.9 


15.7 


70.6 


71.7 


73.4 


Bs 


— ► TjTj 


VI 


6.97 


8.37 


11.6 


10.9 


13.0 


17.7 


56.9 


55.3 


52.8 


B s 




IV 


11.4 


13.2 


17.3 


17.6 


20.4 


26.4 


66.2 


65.6 


64.5 


B s 




I 


4.04 


4.56 


5.70 


4.04 


4.56 


5.70 


0.0 


0.0 


0.0 


B s 




I 


14.8 


16.7 


20.8 


14.9 


16.8 


21.0 


0.9 


0.9 


0.9 


B s 




II 


0.10 


0.003 


0.29 


0.10 


0.002 


0.29 


-1.7 


-36.3 


0.1 


Bs 


-> K°p° 


II 


0.35 


0.04 


1.11 


0.37 


0.07 


1.17 


6.8 


93.8 


5.3 


B s 


-> K°u; 


II,VI 


1.14 


0.16 


1.81 


1.42 


0.26 


1.83 


24.7 


56.7 


1.2 


B s 


-> ^°*r/ 


II,VI 


0.16 


0.13 


0.44 


0.22 


0.21 


0.55 


38.4 


58.2 


24.9 


B s 


-> X°V 


II,VI 


0.08 


0.02 


0.16 


0.10 


0.05 


0.20 


33.6 


131 


21.7 


B s 


-> K+K~* 


IV 


3.05 


3.39 


4.12 


5.03 


5.61 


6.86 


64.7 


65.3 


66.4 


B s 


-> K + *K~ 


IV 


0.89 


0.97 


1.15 


0.90 


0.99 


1.18 


2.2 


2.3 


2.5 


B s 


-► pr? 


V 


0.08 


0.11 


0.25 


0.08 


0.12 


0.25 


1.0 


1.0 


0.8 


B s 


-► pr? 


V 


0.08 


0.11 


0.24 


0.08 


0.11 


0.24 


1.0 


1.0 


0.8 


B s 


— > wry 


V 


0.85 


0.01 


2.60 


1.29 


0.01 


4.15 


51.5 


-1.4 


59.9 


B s 


t 

— > LOT] 


V 


0.84 


0.01 


2.56 


1.28 


0.01 


4.09 


51.5 


-1.4 


59.9 


B s 


1 

-»■ 7T> 


V 


0.13 


0.17 


0.32 


0.13 


0.17 


0.32 


1.9 


1.8 


1.3 


B s 


-> ^ 


VI 


1.85 


0.76 


0.07 


3.78 


1.69 


0.03 


104 


122 


-53.5 


B s 


— ► <f>77 


VI 


0.70 


0.20 


1.49 


1.82 


0.40 


1.14 


161 


107 


-23.5 


B s 


— »■ K^K"* 


IV 


3.24 


4.11 


6.17 


5.52 


6.85 


9.93 


70.6 


66.7 


61.0 


B s 




IV 


0.39 


0.31 


0.18 


0.40 


0.32 


0.19 


0.8 


0.9 


1.0 


Bs 




VI 


0.001 


0.03 


0.30 


0.004 


0.03 


0.40 


118 


1.1 


38.4 


B s 




I 


12.5 


14.1 


17.5 


12.6 


14.2 


17.7 


0.9 


0.9 


0.9 


B s 




II 


0.29 


0.03 


0.94 


0.31 


0.06 


0.99 


6.8 


93.8 


5.3 


B s 




II,VI 


0.24 


0.03 


0.38 


0.30 


0.05 


0.39 


24.7 


56.7 


1.2 


Bs 


-► K+*K~* 


IV 


2.72 


3.02 


3.68 


4.48 


5.00 


6.12 


64.7 


65.3 


66.4 


Bs 




V 


0.15 


0.21 


0.45 


0.15 


0.21 


0.46 


1.0 


0.99 


0.8 


B s 


— ► u;^ 


V 


0.79 


0.01 


2.41 


1.20 


0.01 


3.85 


51.3 


-1.35 


59.9 


B s 




IV 


2.14 


2.71 


4.07 


3.65 


4.53 


6.56 


70.7 


66.8 


61.1 


B s 




VI 


0.03 


0.12 


0.48 


0.05 


0.19 


0.74 


68.5 


58.9 


54.1 


B s 




VI 


17.5 


8.99 


0.42 


29.9 


15.8 


0.98 


71.1 


75.8 


134 
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TABLE IV. B{B S -► PV, VV) (in units of 10~ 6 ) in the SM and model III by using the LCSR 



form factors for B s — > 


or B s 


— > (j) transition, 


and assuming 


k 2 = ml/2, p = 


0.16, 7] 


= 0.34, 


M#+ = 200GeV, 


= 0°, and iV c eff 


= 2, 3, 


oo. 


















SM: B 






Model III: 


B and 5B[%] 




Channel 


Class 


2 


3 


oo 


2 


3 


oo 


2 


3 


oo 


B s -► ^*+vr- 


T 

J. 


4.68 


5.29 


6.61 


4.69 


5.29 


6.61 


0.0 


0.0 


0.0 


5 S -► K+p" 


T 

L 


15.4 


17.4 


21.7 


15.5 


17.5 


21.9 


0.9 


0.9 


0.9 




TT 

_L J_ 


0.12 


0.003 


0.33 


0.11 


0.002 


0.33 


— 1.7 


—36.3 


0.1 


5 S -► 


TT 


0.36 


0.04 


1.16 


0.39 


0.07 


1.22 


6.8 


93.8 


5.3 


5 S -► K°co 


TT VT 

11, V 1 


1.19 


0.17 


1.89 


1.49 


0.27 


1.91 


24.7 


56.7 


1.2 


B s -► if *?? 


TT VT 

11, V 1 


0.18 


0.14 


0.50 


0.24 


0.22 


0.62 


36.7 


57.8 


23.6 


5 S -► if *?/ 


TT VT 


0.09 


0.02 


0.19 


0.11 


0.05 


0.23 


28.7 


137 


17.5 


5 S -► if +if~* 


TV 


3.22 


3.58 


4.35 


5.31 


5.92 


7.22 


64.7 


65.3 


66.4 


5 S -► if +*if~ 


TV 

1 V 


1.04 


1.14 


1.35 


1.06 


1.16 


1.38 


2.2 


2.3 


2.8 


S s -> prj 


V 

V 


0.09 


0.12 


0.26 


0.09 


0.12 


0.26 


1.0 


1.0 


0.8 


Bs OT' 


V 

V 


0.09 


0.12 


0.25 


0.09 


0.12 


0.25 


1.0 


1.0 


0.8 




V 

V 


0.89 


0.01 


2.71 


1.35 


0.01 


4.33 


51.5 


—1.4 


59.9 




v 

V 


0.88 


0.01 


2.67 


1.33 


0.01 


4.27 


51.5 


—1.4 


59.9 


5 S vrV 


v 

V 


0.26 


0.33 


0.63 


0.26 


0.34 


0.64 


1.9 


1.8 


1.3 


E s -> (p7] 


VT 

V J. 


1.36 


0.49 


0.18 


3.04 


1.23 


0.09 


124 


151 


-50.3 


Bs - <M' 


VT 

V 1 


0.38 


0.53 


3.43 


0.87 


0.21 


2.91 


127 


-60.3 


-15.3 


5 S -► if °if°* 


TV 

1 V 


3.42 


4.34 


6.52 


5.83 


7.23 


10.5 


70.6 


66.7 


61.0 


5 S -► if * if 


TV 

± V 


0.46 


0.37 


0.22 


0.46 


0.37 


0.22 


0.8 


0.9 


1.0 


s s if°</> 


VT 

V ± 


0.004 


0.05 


0.36 


0.002 


0.05 


0.50 


-56.3 


1.1 


36.0 


5 S K+*p- 


T 

1 


13.2 


14.9 


18.6 


13.3 


15.0 


18.8 


0.9 


0.9 


0.9 


B s -► if *p 


II 


0.31 


0.03 


0.99 


0.33 


0.06 


1.05 


6.8 


93.8 


5.3 


5 S -► if°*w 


II,VI 


0.26 


0.04 


0.40 


0.32 


0.06 


0.41 


24.7 


56.7 


1.2 


5 S -► if +*if~* 


IV 


2.82 


3.13 


3.79 


4.64 


5.17 


6.33 


64.7 


65.3 


66.4 


B s - pV 


V 


0.27 


0.38 


0.82 


0.28 


0.38 


0.82 


1.0 


1.0 


0.8 


2? s -»■ 


V 


1.43 


0.01 


4.33 


2.16 


0.01 


6.93 


51.5 


-1.4 


59.9 


5 S -► if * if * 


IV 


2.20 


2.80 


4.20 


3.76 


4.67 


6.77 


70.7 


66.8 


61.1 


5 S -► if°*</> 


VI 


0.07 


0.20 


0.66 


0.12 


0.32 


1.03 


68.9 


60.5 


55.1 


5 S # 


VI 


29.9 


15.4 


0.72 


51.1 


27.0 


1.68 


71.1 


75.8 


134 



23 



TABLE V. B(B S -► h x h 2 ) (in units of 10~ 6 ) in model I, with k 2 = ml/2, p = 0.16, rj = 0.34, 
M H+ = 200GeV, tan/3 = 2, and Nf = 2, 3, oo. 



SM: B Model I: B and 



Channel 


Class 


2 


3 


00 


2 


3 


00 


2 


3 


00 


B s -► i^+TT - 


I 


6.33 


7.13 


8.88 


6.34 


7.14 


8.90 


0.1 


0.1 


0.1 




II 


0.19 


0.08 


0.56 


0.19 


0.08 


0.56 


-0.1 


0.0 


0.1 




VI 


0.34 


0.31 


0.78 


0.34 


0.31 


0.79 


0.4 


0.7 


0.4 


£ s -»■ if V 


VI 


0.57 


0.51 


0.76 


0.58 


0.52 


0.77 


1.2 


1.5 


1.3 


B s -»■ K+K~ 


IV 


10.6 


11.7 


14.0 


10.8 


11.9 


14.1 


1.4 


1.3 


1.3 


B s -► vr ?/ 


V 


0.04 


0.06 


0.11 


0.05 


0.06 


0.11 


10.5 


10.0 


7.3 


5 S -► IT T] 


V 


0.04 


0.05 


0.10 


0.05 


0.06 


0.11 


10.5 


10.0 


7.3 


B s -> 7777' 


VI 


13.7 


15.8 


20.5 


13.9 


16.1 


20.8 


1.1 


1.2 


1.3 


5 S -► rj'rj' 


VI 


6.77 


7.48 


9.05 


6.89 


7.63 


9.22 


1.5 


1.5 


1.6 


-B s -> 7777 


VI 


6.95 


8.35 


11.5 


7.03 


8.44 


11.7 


0.8 


0.9 


1.0 


-B s -> 


IV 


11.4 


13.2 


17.2 


11.6 


13.4 


17.5 


1.8 


1.9 


2.0 




I 


4.04 


4.56 


5.70 


4.04 


4.56 


5.70 


0.0 


0.0 


0.0 




I 


14.7 


16.6 


20.8 


14.8 


16.7 


20.8 


0.0 


0.0 


0.0 


-B s -> i^ *7T° 


II 


0.10 


0.003 


0.29 


0.10 


0.003 


0.29 


0.3 


5.1 


0.0 


5 S -► if V 


II 


0.35 


0.04 


1.11 


0.35 


0.04 


1.11 


0.0 


0.5 


0.0 


B s -► if°u; 


II, VI 


1.14 


0.16 


1.81 


1.15 


0.17 


1.81 


0.6 


1.3 


0.0 


-B s -»■ ^°*r/ 


II,VI 


0.16 


0.13 


0.44 


0.16 


0.13 


0.45 


0.4 


0.8 


0.5 




II,VI 


0.08 


0.02 


0.16 


0.08 


0.02 


0.17 


0.6 


2.8 


0.5 


5 S -> K + K~* 


IV 


3.04 


3.38 


4.11 


3.12 


3.45 


4.18 


2.1 


1.9 


1.5 


B s -> K + *K- 


IV 


0.89 


0.97 


1.15 


0.87 


0.96 


1.14 


-1.5 


-1.0 


-0.2 


B s -> /977 


V 


0.08 


0.11 


0.25 


0.09 


0.12 


0.26 


5.4 


5.6 


4.6 


-B s -»■ 


V 


0.08 


0.11 


0.24 


0.09 


0.12 


0.25 


5.4 


5.6 


4.6 




V 


0.85 


0.01 


2.59 


0.86 


0.01 


2.65 


0.8 


-7.4 


2.0 


T"l / 


V 


0.84 


0.01 


2.55 


0.85 


0.01 


2.61 


0.8 


-7.4 


2.0 


5 S -► 7rV 


V 


0.13 


0.17 


0.32 


0.14 


0.18 


0.34 


10.5 


10.0 


7.3 


-Bs -> ^77 


VI 


1.84 


0.75 


0.07 


1.86 


0.76 


0.07 


0.4 


0.2 


2.6 


-Bs -»• 077' 


VI 


0.69 


0.20 


1.49 


0.71 


0.20 


1.51 


1.0 


0.0 


1.2 


5 S -> K°K°* 


IV 


3.22 


4.09 


6.15 


3.30 


4.19 


6.31 


2.0 


2.1 


2.3 


B s -> if°*K° 


IV 


0.39 


0.31 


0.18 


0.39 


0.31 


0.18 


-0.4 


-1.0 


-2.9 




VI 


0.001 


0.03 


0.30 


0.002 


0.03 


0.30 


0.3 


2.1 


1.7 


-B s 


I 


12.4 


14.1 


17.5 


12.5 


14.1 


17.5 


0.0 


0.0 


0.0 


-B s -► if°V 


II 


0.29 


0.03 


0.94 


0.29 


0.03 


0.94 


0.0 


0.5 


0.0 


5 S -► K°*LO 


II, VI 


0.24 


0.03 


0.38 


0.24 


0.04 


0.38 


0.6 


1.3 


0.0 


B s -► K+*K-* 


IV 


2.71 


3.02 


3.66 


2.78 


3.08 


3.73 


2.1 


1.9 


1.5 


B s - pV 


V 


0.15 


0.21 


0.45 


0.16 


0.22 


0.47 


5.4 


5.6 


4.6 


S s -> u;^ 


V 


0.79 


0.01 


2.40 


0.80 


0.01 


2.46 


0.7 


-7.4 


2.0 


5 S -► if°*K°* 


IV 


2.13 


2.70 


4.06 


2.17 


2.77 


4.17 


2.0 


2.1 


2.3 


-B s -► if°*0 


VI 


0.03 


0.12 


0.48 


0.03 


0.12 


0.49 


3.8 


2.9 


2.4 


5s # 


VI 


17.4 


8.95 


0.42 


17.7 


9.08 


0.42 


1.1 


1.1 


0.5 
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TABLE VI. B(B S -► h x h 2 ) (in units of 10~ 6 ) in model II, with k 2 = m 2 b /2, p = 0.16, r? = 0.34, 
M H + = 200GeV, tan/? = 2, and iV c cfr = 2, 3, oo. 



SM: B Model II: B and &B[%] 



Channel 


Class 


2 


3 


00 


2 


3 


00 


2 


3 


00 


B s -> iT+TT - 


I 


6.33 


7.13 


8.88 


6.29 


7.09 


8.84 


-0.6 


-0.6 


-0.6 


-B s -> if°7r° 


II 


0.19 


0.08 


0.56 


0.18 


0.07 


0.55 


-4.7 


-13.4 


-2.9 


5 S -► K u r] 


VI 


0.34 


0.31 


0.78 


0.31 


0.28 


0.74 


-7.5 


-9.8 


-5.2 


B s -»■ if ?/ 


VI 


0.57 


0.51 


0.76 


0.52 


0.45 


0.70 


-9.5 


-11.6 


-9.1 


£ s -»■ K+K~ 


IV 


10.6 


11.7 


14.0 


9.58 


10.5 


12.5 


-10.2 


-10.4 


-10.8 


B s -► 7r°r/ 


V 


0.04 


0.06 


0.11 


0.05 


0.06 


0.11 


10.5 


10.0 


7.3 


5 S -► IT T] 


V 


0.04 


0.05 


0.10 


0.05 


0.06 


0.11 


10.5 


10.0 


7.3 


B s -> 7777' 


VI 


13.7 


15.8 


20.5 


12.1 


14.0 


18.2 


-11.7 


-11.6 


-11.4 


5 S -► 7/7/ 


VI 


6.77 


7.48 


9.05 


5.94 


6.56 


7.91 


-12.5 


-12.7 


-12.9 


-B s -> W 


VI 


6.95 


8.35 


11.5 


6.21 


7.49 


10.4 


-10.9 


-10.6 


-9.9 


-B s -> 


IV 


11.4 


13.2 


17.2 


10.3 


11.9 


15.6 


-11.9 


-11.7 


-11.4 




I 


4.04 


4.56 


5.70 


4.04 


4.56 


5.70 


0.0 


0.0 


0.0 


5 S -► K+p- 


I 


14.7 


16.6 


20.8 


14.7 


16.6 


20.8 


-0.2 


-0.2 


-0.2 


£ s -> K°*ir° 


II 


0.10 


0.003 


0.29 


0.10 


0.003 


0.29 


0.8 


16.5 


0.0 


B s -► if °p° 


II 


0.35 


0.04 


1.11 


0.35 


0.03 


1.10 


-1.1 


-16.7 


-1.0 


5 S -► if w 


II,VI 


1.14 


0.16 


1.81 


1.09 


0.15 


1.80 


-4.5 


-10.3 


-0.2 


-B s -»■ ^°*r/ 


II,VI 


0.16 


0.13 


0.44 


0.15 


0.12 


0.42 


-7.3 


-11.1 


-4.7 




II,VI 


0.08 


0.02 


0.16 


0.07 


0.02 


0.16 


-4.7 


-17.9 


-2.9 


5 S -> K + K~* 


IV 


3.04 


3.38 


4.11 


2.74 


3.03 


3.66 


-10.3 


-10.7 


-11.2 


B s -> K + *K- 


IV 


0.89 


0.97 


1.15 


0.87 


0.95 


1.14 


-2.0 


-1.6 


-0.8 


B s -> ,977 


V 


0.08 


0.11 


0.25 


0.09 


0.12 


0.26 


5.4 


5.6 


4.6 


-B s -»■ W 


V 


0.08 


0.11 


0.24 


0.09 


0.12 


0.25 


5.4 


5.6 


4.6 


£ s -> W77 


V 


0.85 


0.01 


2.59 


0.77 


0.01 


2.33 


-9.9 


-7.4 


-10.3 


T"l / 


V 


0.84 


0.01 


2.55 


0.76 


0.01 


2.30 


-9.9 


-7.4 


-10.3 


5 S -► 7rV 


V 


0.13 


0.17 


0.32 


0.14 


0.18 


0.34 


10.5 


10.0 


7.3 


-Bs -> 0r? 


VI 


1.84 


0.75 


0.07 


1.51 


0.60 


0.09 


-18.3 


-20.5 


26.8 


-Bs -»• <fnf 


VI 


0.69 


0.20 


1.49 


0.55 


0.21 


1.60 


-21.3 


7.5 


7.6 


B s -> K°K°* 


IV 


3.22 


4.09 


6.15 


2.85 


3.64 


5.54 


-12.0 


-11.3 


-10.1 


B s -> K°*K° 


IV 


0.39 


0.31 


0.18 


0.39 


0.31 


0.18 


-0.6 


-1.3 


-3.2 


B s 


VI 


0.001 


0.03 


0.30 


0.002 


0.03 


0.28 


6.6 


2.0 


-6.6 


-B s if +*/>- 


I 


12.4 


14.1 


17.5 


12.4 


14.0 


17.5 


-0.2 


-0.2 


-0.2 


B s -► if °V° 


II 


0.29 


0.03 


0.94 


0.29 


0.03 


0.93 


-1.1 


-16.7 


-1.0 


5 S -► if °*LJ 


II,VI 


0.24 


0.03 


0.38 


0.23 


0.03 


0.38 


-4.5 


-10.3 


-0.2 


-B s -► if +*if~* 


IV 


2.71 


3.02 


3.66 


2.43 


2.70 


3.27 


-10.3 


-10.7 


-11.2 


B s - pV 


V 


0.15 


0.21 


0.45 


0.16 


0.22 


0.47 


5.4 


5.6 


4.6 


S s -> u;^ 


V 


0.79 


0.01 


2.40 


0.71 


0.01 


2.16 


-10.0 


-7.4 


-10.3 


5 S -► if * if * 


IV 


2.13 


2.70 


4.06 


1.88 


2.41 


3.66 


-12.1 


-11.3 


-10.2 


-B s -► if°*</> 


VI 


0.03 


0.12 


0.48 


0.03 


0.11 


0.44 


-10.0 


-9.2 


-8.8 


5s # 


VI 


17.4 


8.95 


0.42 


15.2 


7.75 


0.33 


-13.0 


-13.8 


-21.4 
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TABLE VII. CP-violating asymmetries Acp(B s — > h\h,2) (in percent) in the SM and model 
III, with k 2 = ml/2, p = 0.16, rj = 0.34, M H + = 200GeV, 6 = 0°, 30°, and Nf = 2, 3, oo. 









SM 




Model III: 9 


= o°, 


Model III: = 


30°, 


Channel 


CP-Class 


2 


3 


00 


O 

z 


Q 
O 


OO 


O 

z 


9 



00 


B s -> i^ + 7T~ 


1 


10.2 


10.2 


10.3 


9.95 


n nc 

9.9o 


9.97 


1 n 
1U.Z 


1 n 

10.2 


10.3 




2 


-1.99 


-3.98 


4.50 


—3.18 


—3.87 


4.48 


r'C\ 

—2.60 


A El 

—4.51 


4.36 


s SI 


2 


-4.73 


-3.62 


2.93 


A Q/l 


9 E.Q 

— o.oy 


2.17 


A 7fi. 
— 4. ( D 


/I in 

— 4. iy 


1.67 


B s -► 


2 


0.31 


-2.86 


-4.74 


n En 
— U.oU 


nn 

— z.yy 


-4.96 


1 A A 

— 1.44 


9 Q/l 

— o.o4 


-4.77 


5 S -► K + K~ 


2 


-1.71 


-1.74 


-1.77 


1 /in 
— 1.4U 


1 /i 1 

— 1.41 


-1.43 


O 9 Q 
— Z.OO 


a n 
— Z.4U 


-2.45 


B s -> vr°r/ 


2 


-2.72 


-0.24 


2.93 


fin 
—z.oy 


n o/i 
— U.Z4 


2.92 


fin 

— z.oy 


n O/i 
— U.Z4 


2.92 


5 S -► vrV 


2 


-2.72 


-0.24 


2.93 


— z.oy 


n o/i 

— U.Z4 


2.92 


O fiQ 

— z.oy 


n o/i 

— U.Z4 


2.92 


-B s — > 7777' 


2 


0.07 


0.05 


0.01 


n nE 
U.Uo 


n HQ 
U.Uo 


0.01 


1 no 
— l.Uo 


1 nn 

— i.uy 


-1.11 




2 


-0.16 


0.03 


0.36 


n 1 9 
— U.lo 


n no 
U.Uz 


0.27 


1 9 /I 

— l.o4 


1 01 

— l.zl 


-0.98 


— ► Tiri 


2 


0.30 


0.06 


-0.31 


U.Z4 


n nE 
U.Uo 


-0.25 


n qi 
— U.ol 


n n7 

— u.y / 


-1.24 




2 


0.05 


0.04 


0.04 


U.Uo 


n hq 
U.Uo 


0.03 


1 1 n 

— i.iy 


1 1 Q 

— l.lo 


-1.17 


5 s _> K* + TT- 


1 


0.34 


0.34 


0.34 


n 9 A 
0.34 


n /i 
U.34 


0.34 


n 9 /i 
U.34 


n 9 /i 
U.o4 


0.34 


B s -> 


1 


5.56 


5.56 


5.55 


e eo 
5.52 


o.ol 


5.51 


5.6U 


5.0U 


5.59 




1 


-8.45 


-26.6 


6.17 


—8.59 


/i 1 
—41.2 


6.17 


q An 
— 8.40 


A n 1 
— 4U.1 


6.27 


5 S if Sp° 


1 


-22.6 


-16.3 


15.5 


—21.1 


T nn 

—7.99 


14.8 


on 9 
— 20.O 


— 16.6 


15.4 


B q -> KSu; 


2 


3.88 


-1.60 


4.53 


9 1 e 


— z.Ul 


4.65 


z.7o 


nn 

— z.yU 


4.64 




1 


-27.3 


-0.36 


23.1 


— ly.o 


— U.U1 


18.6 


01 a 
— zl.o 


— 4.DU 


18.7 




1 


52.8 


30.8 


-45.0 


on n 
oy.y 


14.1 


-36.7 


4z.o 


nn 

— z.uy 


-35.7 


5 S -> K+K-* 


1 


-30.5 


-30.9 


-31.7 


— lo.o 


inn 

— ly.u 


-19.4 


oo n 
— zo.y 


/i 

— z4.z 


-24.6 


B s -> K + *K~ 


1 


1.40 


1.45 


1.53 


1.37 


1 An 


1.49 


1.37 


1 /i 1 
1.41 


1.49 


B q — > or? 


2 


-3.03 


-0.26 


3.02 


no 

—3.02 


n or* 

— 0.26 


3.01 


no 

— 3.U2 


n or* 

— 0.26 


3.01 


5, -»■ or?' 


2 


-3.03 


-0.26 


3.02 


9 09 
— O.UZ 


n ofi 

— U.ZD 


3.01 


9 no 
— o.uz 


n ofi 
— u.zo 


3.01 


£? s — > CJ?7 


2 


-0.87 


-1.01 


-0.84 


n 71 
— U. / 1 


1 no 
— l.Uz 


-0.68 


1 fi7 
— 1.0/ 


1 no 
— l.Uz 


-1.74 


S s — > Wt/ 


2 


-0.87 


-1.01 


-0.84 


n 71 
— U. ( 1 


1 no 
— l.UZ 


-0.68 


1 fi7 
— 1.0/ 


1 no 
— l.UZ 


-1.74 


5 S -► 7T°0 


2 


-2.72 


-0.24 


2.93 


fin 
— z.by 


n o/i 
— U.z4 


2.92 


fin 

—z.oy 


n o/i 
— U.Z4 


2.92 


6, — ► (irj 

a i t 


2 


0.49 


0.17 


2.87 


n 71 
— 0. / 1 


1 no 
— 1.02 


-0.68 


1 fi7 

— 1.6/ 


1 no 
— 1.02 


-1.74 


B s -> 


2 


-0.31 


0.25 


-0.74 


n 71 
— U. / 1 


1 no 
— l.Uz 


-0.68 


1 fi7 
— 1.0/ 


1 no 
— l.Uz 


-1.74 




1 


-1.25 


-1.21 


-1.13 


1 no 
— l.Uz 


n no 

— u.ys 


-0.93 


fi no 
— o.yo 


P. A P. 
— 0.40 


-5.72 




1 


-0.02 


-0.03 


-0.04 


n no 
— 0.02 


n no 
— 0.03 


-0.04 


n no 
— U.U2 


n no 
— 0.03 


-0.04 




2 


-3.38 


-3.45 


-3.27 


O QQ 
— Z.OO 


Q /I E 

— o.4o 


-3.30 


A A O 

— 4.4Z 


9/17 

— 0.4/ 


-3.84 




1 


5.56 


5.56 


5.55 


r r 

5.52 


5.51 


5.51 


5.60 


5.60 


5.59 




1 


-22.6 


-16.3 


15.5 


-21.1 


-7.99 


14.8 


-20.3 


-16.6 


15.4 




1 


25.5 


13.1 


5.43 


20.6 


8.57 


5.37 


20.9 


4.39 


5.46 


5 S -► K+*K-* 


3 


-3.86 


-3.87 


-3.90 


-3.25 


-3.27 


-3.29 


-4.05 


-4.07 


-4.10 


Bs - pV 


3 


-3.03 


-0.26 


3.02 


-3.02 


-0.26 


3.01 


-3.02 


-0.26 


3.01 


£ s -> u;^ 


3 


-0.87 


-1.01 


-0.84 


-0.71 


-1.02 


-0.68 


-1.67 


-1.02 


-1.74 


5 S -► i^°*i? * 


3 


0.05 


0.05 


0.04 


0.04 


0.03 


0.03 


-1.18 


-1.13 


-1.06 


5 S -► 


1 


4.41 


3.46 


3.01 


2.90 


2.40 


2.15 


-3.03 


-2.27 


-1.91 


5, -► # 


3 


0.05 


0.06 


0.12 


0.04 


0.04 


0.06 


-1.19 


-1.24 


-1.79 
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TABLE VIII. 






) (in percent) in 


models I and II. 


with k 2 = ml/2, p = 


- U.1D, 


V = 


= 0.34, M H+ = 


- ZUUvjc V , 


ran /j — 


2, and Nf = 2, 


3, 00. 




















SM 






Model I 




Model II 




Channel 


CP-class 


2 


3 


00 


2 


3 


00 


2 


3 


OO 


B a 


-► K + tt- 


1 


10.2 


10.2 


10.3 


10.2 


10.3 


10.3 


10.3 


10.3 


10.3 


Bs 




2 


-1.99 


-3.98 


4.50 


-1.98 


-3.98 


4.79 


-1.63 


-4.01 


4.84 


B s 


-► K° s r] 


2 


-4.73 


-3.62 


2.93 


-4.74 


-3.62 


2.92 


-4.65 


-3.63 


3.11 


B s 




2 


0.31 


-2.86 


-4.74 


0.28 


-2.86 - 


-4.75 


0.52 


-2.82 


-4.64 


B s 


-► K+K- 


2 


-1.71 


-1.74 


-1.77 


-1.70 


-1.72 - 


-1.76 


-1.80 


-1.82 


-1.86 


Bs 


— > TT°rj 


2 


-2.72 


-0.24 


2.93 


-2.59 


-0.23 


2.84 


-2.59 


-0.23 


2.84 


B s 


-virV 


2 


-2.72 


-0.24 


2.93 


-2.59 


-0.23 


2.84 


-2.59 


-0.23 


2.84 


B s 


— > 7777' 


2 


0.07 


0.05 


0.01 


0.07 
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FIG. 1. Branching ratios B(B S — > (fir]) versus M H + and l/Nf in the SM and models II 
and III by using the BSW form factors. For (a) and (b), we set Nf = 3 and M H + = 200GeV, 
respectively. The four curves correspond to the theoretical predictions in the SM (dotted line), 
model II (dot-dashed curve), model III with 6 = 0° (solid curve) and 6 = 30° (short-dashed 
curve), respectively. 
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FIG. 2. Same as Fig.| but for the decay B s -> K + K~ . 
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FIG. 4. Same as FigJT] but for the decay B s -> K+R-*. 
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FIG. 6. CP-violating asymmetries Acp of B s — ► K + K decay versus and 1/N° 

in the SM and models II and III. For (a) and (b), we set Nf = 3 and M H + = 200GeV, 
respectively. The four curves correspond to the theoretical predictions in the SM (dotted line), 
model II (dot-dashed curve), model III with 6 = 0° (solid curve) and 6 = 30° (short-dashed 
curve), respectively. 
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FIG. 7. Same as Fig.| but for decay B s -> K + K~*. 
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FIG. 8. Branching ratios and CP-violating asymmetries of B s — ► K°*r)' decay versus 
l/iV c cfr in the SM and models II and III, assuming M H + = 200GeV and tan/3 = 2. The 
four curves correspond to the theoretical predictions in the SM (dotted curve ), model II 
(short-dashed curve), model III with 9 = 0° (long-dashed curve) and 9 = 30° (solid curve), 
respectively. 
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